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Abstract 

Using a recursive method we construct dimer and nondimer variational ansatzs of the 
ground state for the two-legged ladder, and compute the number of dimer coverings, 
the energy density and the spin correlation functions. The number of dimer coverings 
are given by the Fibonacci numbers for the dimer-RVB state and their generalization 
for the nondimer ones. Our method relies on the recurrent relations satisfied by the 
overlaps of the states with different lengths, which can be solved using generating 
functions. The recurrent relation method is applicable to other short range systems. 
Based on our results we make a conjecture about the bond amplitudes of the 2- 
ladder. 

PACS numbers: 05.50.+q, 75.10.-b, 75.10.Jm 
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I. INTRODUCTION 



There is an increasing theoretical and experimental interest in systems formed by a finite 
number ni of coupled chains, currently known as ladders (for a review see Theoretical 
studies suggested that antiferromagnetic spin ladders should be gapped (gapless) depending 
on whether rii is an even (odd) number [0], 0, [^], 0, This prediction has been 

confirmed experimentally in compounds like (VO)2P207 [rii = 2) 0, although some doubts 
has been recently casted on the ladder structure on this material |TD[ , SrCu203 {rii = 2) [|TT 



Sr2Cu305 {rii = 3), etc. Moreover, the theory also predicts that upon doping the even-n/ 
ladders should become superconductors due to the pairing of holes 0,11], [|12| . The recent 
discovery of superconductivity in the 2-legged ladder compound Sro.4Cai3.6Cu2404i under 
high pressure |jl3| may perhaps constitute a confirmation of this prediction [Q. On the 
other hand the odd-ri; ladders are not expected to superconduct. 

Altogether the ladders fall in two different universality classes depending on the 
even/oddness of ni, which in the limit ni ^ oo should converge to the same class. An 
example of this is given by the behaviour of the spin gap of the even-rz; ladders which 



vanishes exponentially with ni |T5[, [T^, ||T7|. 

A field theoretical characterization of these two universal classes can be obtained by 
mapping the spin ladders into the 1 + 1 non- linear sigma model |]TB|, ITB]. The value of 
the coupling constant 6, which multiplies the instanton number in the action, is given by 



6 = 2tx Srii, where S is the spin of the chain |jT9|], JTSl, [^. For S* = 1/2 and rii even one gets 
^ = (mod 27r), which corresponds to a sigma model with a dynamically generated gap, 
while for rii odd one gets 6 = tt (mod 27r), which corresponds to a gapless sigma model which 
fiows under the RG to the SU{2) level 1 Wess Zumino model [0. Thus the behaviour of 
ladders parallels that of the spin chains as function of the spin, as first cojectured by Haldane 
using precisely the mapping of the spin chains into the non-linear-sigma model p2[, ||23|| . 



There is an alternative explanation of the qualitative difference between the even/odd 
ladders based on the RVB (Resonating Valence Bond) theory. According to the authors of 
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ref . , the even-n/ ladders are short ranged RVB systems with confinement of spin deflects, 
which leads to the existence of a spin gap and exponential decaying correlation functions; 
while the odd-ni ladders are long-range RVB systems with no confinement of deflects, no 
gap and power-law correlations. Strictely speaking, the above RVB interpretation based 
on the original ideas of Liang et al. |^ concerning the 2D AF-Heisenberg system, are an 



intuitive picture to explain the numerical results obtained using the DMRG |25| (Density 
Matrix Renormalization Group). It is therefore interesting to test the RVB picture using 
diflerent techniques in order to confirm and get further insights into the short range nature 
of the even legged-ladders. 

The first step in this direction is the study of the dimer-RVB state (also termed NNRVB, 
standing for nearest-neighbour RVB). This has already been pursued in ||26|, for 
ni = 2. Using the DMRG, the main conclusion of is that one has to consider valence 
bonds between sites which are not NN (nearest-neighbours), going in that way beyond the 
dimer-RVB towards a short-range RVB state, whose structure has not yet been studied in 
detail. 

From a mathematical viewpoint the dimer-RVB states are relatively easy to handle due to 
the work of Sutherland which gives an elegant diagrammatic way of computing overlaps 



between the dimer states which form the dimer-RVB state. However the combinatorics 
of the nondimer-RVB states has not yet been worked out and their construction seems a 
priori a difficult task. We shall show in this paper that there is a way to circumvent the 
computational problem associated with short-range RVB states, which is based on the use 
of recurrent relations. The main idea is to build up the ground state of a ladder with N + u 
rungs using the knowledge of the ground states with N, N + 1, . . . , N + u — 1 rungs. This 
is achieved by a recurrent relation (RR) which gives the ground state \N + u) in terms of 
the ground states \N + u — 1) , . . . , \N). We shall call u the order of the recurrent relation. 
The "matching" of the various g.s.'s within the RR is achieved by means of a collection 
of z/ states . . . \(pu), which are the elementary building blocks of our method. Loosely 
speaking, is a state which contains at least a bond of length k. For even-n/ ladders. 



with ^ a few lattice spacings, we expect an adequate description of the ground state for 
small values of v. The use of RR's allow also the introduction, in a natural fashion, of 
variational parameters which we can determine by the standard minimization procedure of 
the ground state energy. The whole approach based on RR's is analytic, powerful and rather 
straightforward, and we believe it represents a significant methodological improvement as 
compared with previous analytic approaches to the study of RVB states, which is applicable 
to other type of short range systems. 

The organization of the paper is as follows. In section II we introduce the recurrent 
relations and compare the states generated by them with the conventional RVB states. In 
section III we deduce the RR's for the norm of the ground state and the expectation value 
of the Hamiltonian. These RR's are solved in section IV using generating functions, which 
allow us to obtain quite easily an analytic formula for the ground state energy density. The 
results of the minimization of the ground state energy respect to the variational parameters 
of the ansatzs corresponding to v — 2 and 3, are presented in section V. In section VI we 
derive analytic expressions for the correlation length ^ and give their numerical values for 
V — 2 and 3. Finally in section VII we discuss our results and the perspectives of the RR 
method. 
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II. RVB AND RVA STATES 



Generically, an RVB state is a linear superposition of singlets states constructed by 
pairing up the N spins of a system into y bonds (A^ is even). Making an analogy with 
the BCS states one can think of an amplitude h{i — j) of creating a bond, denoted by (zj), 
between the sites i and j p9[. We shall work with bipartite lattices and such that i (j) 
belongs to the sublattice A (B). The total RVB wave funtion is given by p^ , 



\i')RVB = ^(^1 ~ ^l) ■ ■ ■ - jN/2) (iljl) ■ ■ ■ {iN/2iN/2) (l) 

ieA.'jeB 

where the amplitudes h{i—j) are all positive or zero in order to satisfy the Marshall sign rule 
30| . For a dimer state h{i — j) 7^ unless i and j are nearest neighbours (NN). In a short- 



range RVB state, h{i — j) decreases exponentially with the distance while in a long-range 
RVB state, h{i — j) decays algebraically as \i — for some power p. 

Let us consider the RVB state (|lD for the two- leg ladder depicted in figure 1. We have 
labeled the sublattice A hj k = 1,2, . . . , N and the sublattice B by A; = 1, 2, . . . , A^. The 



general ansatz (|l|) only depends on the amplitudes hi of the bond {k, k + I) for / = 0, 1, . . .. 
Let us consider some examples in order to get familiar with the notation, while going through 
the details of the formalism. 

A. Columnar State 

Let us first assume that only /iq is non zero. We normalize it as /iq = 1- In this case the 
corresponding RVB state for a ladder with rungs takes the simple the form, 

N _ 

\N) = \{{k,k) (2) 



k=l 



This state is actually the exact ground state of the rung Hamiltonian defined by, 

N 

Hrnng = J Sfc ■ S-^ (3) 

k=l 

where the vertical coupling constant J' is antiferromagnetic (i.e. J' > ). 
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FIG. 1. A two-legged ladder showing the labels corresponding to sublattices A and B. 

The correlation length C, of the state measured by the spin-spin correlator along the 
legs is exactly 0. 

A trivial observation is that the states + 1) and |A^) given by eq. (0) are related by the 
equation, 

|iV+l) = |iV)®|0i)^+i (4) 

where |0i)Ar+i denotes the singlet state located at the {N + 1)— rung, namely, 

|0i)jv+i = (iV + l,iVn) (5) 

Eq.(||) is a 1** order RR (i.e. z/ = 1), which upon — 1 iterations reproduces the state 
given the initial condition |1) = \4>i). 



B. Dimer State 

Let us assume that = 1 and hi are both different from zero. The variational parameter 
hi gives the amplitude of a horizontal bond in the RVB ansatz (|1|), i.e. 



\N) = Y.h-,^il,ki)i2,k,)---iN, k 



N) 



(6) 



where ni is the number of horizontal bonds. The state is a linear superposition of dimer 
states of the form depicted in figure 2, and for this reason is called a dimer-RVB state. 



Observe that all the horizontal bonds come in pairs, for each horizontal bond at one leg 
forces the presence of a companion in front of it at the opposite leg. 



FIG. 2. A typical dimer state in a 2-legged ladder. 



It is apparent that upon switching the Hamiltonian which contains the couplings along 
the legs of the ladder, namely 

Hleg = J'^iSk- Sj^ + Sfc+i • S^) (7) 

fc=l 

the dimer states (§) will become more probable than the columnar state @ which contains 
only vertical rungs, and for J = J' one expects the amplitude hi to be close to 1. This 
expectation is based on the resonant mechanism that motivates the whole RVB approach 



| 3T| , and which is the main cause of the substantial lowering of the energy for the RVB 
states (see figure 3). 

The key observation for our purposes is that the dimer state @ of a ladder with open 
boundary conditions can be generated from a 2"''^ order RR given by: 

\N + 2) = \N+1)® \4>i)n+2 + u\N)® \(f)2)N+l,N+2 (8) 

where the state denoted by \4>2) n+i,n+2 is made up of a pair of horizontal bonds located 
between the rungs at (A^ + 1, + 2), i.e.. 



|02)iv+i,7v+2 = {N + 1,N + 2){N + 2,N + 1) (9) 
Comparing with (§) we obtain the relation between the parameters hi and u, 

u = hl (10) 
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FIG. 3. The basic bond resonance mechanism for an elementary plaquette. 



We refer to figure 4 for a diagrammatic selfexplanatory representation of the RR (|^). 

One expects the dimer state (|^) to describe correctly ground states for which the corre- 
lation length is at most one. 



N+2 



N+1 



I . u [7] 



FIG. 4. Diagrammatic representation of the Recurrent Relation for u = 2. 



Fan and Ma have computed the ground state energy and the number of dimer states 
contained in the state |A^) [^. Let us call that number for reasons that will become 
clear below. F/v is given by, 

1 



do 



(11) 



where ai^2 = |(1 =t V^)- The key fact is to recognize eq.([Tl|) as the Binet's formula for the 
j\j-th pibonacci number ||3^. This last result inmediately follows from the RR satisfied by 
F/v, as a consequence of eq.(||), namely, 



"^+2 — -^AT+l + 



(12) 



This equation together with the initial values. 



Fo = Fi = 1 



(13) 



reproduce the well known Fibonacci sequence. 
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1,1,2,3,5,8,... 



(14) 



Another way to arrive to this resuh is by counting the number of dimer states in (^. 
Calhng M the number of pairs of horizontal bonds ( M = 0, 1, . . . , [N/2], where [x] denotes 
the integer part of x), one gets. 



[7V/2] 

fn=y: 

M=0 



V 



M 



(15) 



This formula is well known in the theory of partitions [^]. 



C. Nondimer States 

These states are obtained whenever hi is non-zero for bonds connecting non NN (nearest- 
neigbour) states. As an example, let us assume that the non vanishing amplitudes are = 1, 
hi and h2. In figure 5 we show "local" structures formed with these bonds. The number of 
configurations of this type increases enormously with the number of rungs. However not all 
of them are independent. The dimer states are all linearly independent. 






a) 



b) 



c) 





d) 



e) 



FIG. 5. A picture of several non-dimer states showing their increasing variety. 



9 



For example, the state depicted in fig. 5c) can be written as a linear combination of dimer 
states and the states in fig. 5a) and fig. 5b) (see figure 6). 




FIG. 6. A non-dimer state as a linear combination of the 5 linearly independent states defining 
the V = 2> Recurrent Relation. 

It is clear from figures 5 that in order to generate the nondimer RVB states in full 
generality one has to consider RR's with an arbitrary high order v. Thus, configurations 
5a), 5b) and 5c) require a 3'"'^ order RR, configuration 5d) a 4*^^ order configuration, and 
so on. For physical reasons one expects /i2 < hQ,hi, so that configurations 5d) and 5e) are 
much less probable than configurations 5a),b),c). Hence, if we restrict ourselves to these 
latter class of configurations we can again generate that class iteratively by means of the 
following RR, 

\N + 3) = \N + 2)® |0l)jV+3 + U\N + 1) ® |02)iV+2,7V+3 + V\N) ® |03)jV+l,iV+2,iV+3 (16) 




FIG. 7. A pictorical representation of the non-dimer state l^s) defining the v = 3 Recurrent 
Relation. 
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where |03)Ar+i,Ar+2,Af+3 is pictured in figure 7 and its precise expression reads, 

\h)i,2,3 = ^[(1,2)(2,3)(3J) + (1,3)(2,T)(3,2)] (17) 

The prefactor -y/l is for normalization. The last term in (|T6|) generates states with two 
horizontal bonds of type hi and one bond of type h2, which follows the chess knight's move 
||35| ( this kind of bonds have been recently considered in the study of the Hubbard model 
on 2 X 2 X 2 and 4x4 lattices in ref. ). The relation between hi, /12 and v is given by, 



^jlv = h\ h, (18) 

We have not included in l^s) the configuration 5c) since it depends linearly on the 
remaining configurations contained in ([161) . In a sense eq.(|16D generates recurrently the 



most important RVB configurations with bonds of type ho, hi and /i2- The RR (|T6D will 
be used to generate all the states |A^) with > 1 provided we make the following formal 
identifications, 

I 1, A^ = 

\N) = (19) 
I 0, N <0 

These are also the initial conditions of the RR (|]). 

The states generated recurrently from equations similar to (|^) and (|T6D shall be called 
hereafter RVA states (standing for Recurrent Variational Approach), to distinguish them 
from the RVB states of the form (|l|). We have proved above that the RVA states produced 
by 1*^* and 2"^^ order RR's actually coincide with RVB states, but this is not true for higher 
order RVA states. The main advantage in working with RVA states is that they can be 
treated analytically. 

It is quite clear that we can perform some generalizations of the previous ideas. The 
general form of a RR of order u can be written as, 

V 

(20) 
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where are normalized states which must be chosen to be hnearly independent from the 
states generated in the previous steps, while ui are variational parameters. 

The RR satisfied by the number, of linearly independent states generated by (pO[), 
is given by: 

V 

Fn+u = ^N+u-l (21) 
1=1 

subject to the initial conditions Fq = 1 and Fjv = for < 0. Therefore the higher order 
RR's correspond to generalizations of the Fibonacci numbers. In section IV we shall get the 



following Binet's formula for 



N 1 



F^N^ = i:\{^^ (22) 

i=i j^i 

where are the roots of the polynomial y" — y"^^ — . . . — y — 1. 



III. RECURRENCE RELATIONS FOR OVERLAPS 

In this and the following section we shall compute the value of the energy 
{N\H\N)/{N\N) of the RVA states \N). 

Let us proceed progresively and first consider the 2""^ order RR given by eq.(P). It is 
convenient to introduce the quantities and Yn as: 

= {N\N) ^^^^ 
Yn= n{<Pi\(^{N-1\N) 

Now one can easily derive the following RR's, 

Zn+2 = Zn^i + uY^^i + u'^Zn 

(24) 

Yn+2 = Zn^i + 
where we have made use of the result, 

n+i{4>i\4>2)n,n+i = ^\4>i)n (25) 



The RR's (^J) together with the initial conditions, 
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determine Z^- and for arbitrary values of A^. This will be done in the next section using 
generating-function methods. 

The Hamiltonian of a ladder with rungs, is given by 

Hn = Hicg + i^rung (27) 

where H\eg and -ffrung are defined in (H) and respectively, The RR method applied to 
Hamiltonian overlaps requires the following definitions. 

En = {N\Hn\N) 

Dn= N{(t>i\®{N ~1\Hn\N) 



A straightforward computation using eqs. (H) and (pSD leads to the following RR's for 
and Dni 

En+2= En+i + J'eoZN+i + u{DN+i + i2J+ J')eoYN+i) + u\EN + 2JeoZN) 

(29) 

Dn+2 = Ejy+i + J'eoZjy^i + ^{D]\r+i + (2 J + J')^o^Af+i) 
where eo = —3/4 is the lowest eigenvalue of the operator 81-82. 
The initial conditions for En and Dn are. 



(30) 



Eq = 0, El = J'tQ 
Do = 0, Di = J'eo 

Before showing the power of the RR method by computing the values of several physical 



quantities, we set off for the generalization of equations (|23D and p8| ) to arbitrary values of 
the order u of the RR. Let us define the following quantities: 

Zn,1 =N,...,N-l+l ® {N-l\N) 

(/ = 0,1,...,J/-1) (31) 

En,i =N,...,N-i+i {4>i\ ® (A^ - 1\Hn\N) 
where we identify Zn,o = {N\N) and En,o = {N\Hn\N). 
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The RR's (pO]) for the ground states |A^), does not automatically imply that the overlaps 
Zjv/ also satisfy RR's. However this is guaranteed provided the states satisfy the 

following eqs. 



N,...,N-k+l\(Pk\(Pl) N-l+l,...,N 



kl \Vl-k) N-l+l,...,N-k 



I > k 



(32) 



Qu = l l = k 

^kl N-l,...,N~k+l{4>k-l\ I < k 

As a matter of fact, eq.(P3[) gives an example of conditions (^). The states |02) and 
103) also satisfy eqs. (0), with the overlapping matrix Vl given by, 

1 I2\ 



^(.=3) 



/ 



2 2 1 
VV5 V5 ^ 



(33) 



Eqs.([33|) can be given a geometrical meaning in terms of the disentangling of connected 
bonds ( see figure 8). 

For z/ = 3, beside the overlaps defined in eqs.(|3TD, we need in addition the following 
matrix element, 



{N -2\N) 



(34) 



The RR's for the overlaps Zj^^i{l = 0, 1, 2) and can be derived from eq.(|TBp, 

Zn+3 = Zn+2 + U^Zn+I + v'^Zn + 2uVLi2Zn+2,1 + '^vVLi-iZN+2,2 + '^UvVL2-iZN+l,l 

Zn+-1,,1 = Zn^2 + uVLi2Zfq^2,l + vVLi-iZN^2,2 

Zn+3,2 = uZjsf^i + ^12Zn^2,1 + V^23Zn+1,1 

Wn+3 = uQuZn+I + Zn+2,1 + ^^^12^13^Af+l,l 

The initial data to solve these RR's are, 

Zq = 1 Zq^i = Zo,2 = 

Zl = 1 ^1,1 = 1 ^1,2 = 



(35) 



Z2 = l+U + U^ ^2,1 = 1 + f ^2,2 = I 



Wo = 

Wi = 

U W2 = l + ^ 



(36) 
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= 1/2 



1/2 



*3> = ^(2/5) ( 




+ 




) 



V(2/5) 



V(2/5) (|). 



02 O = V(2/5) ( 




+ 




) 



V(2/5) 



V(2/5) \ 



FIG. 8. Graphical representation of the overlapping calculus for the building block states. 



In order to derive RR's for the overlaps of the Hamiltonian Hn it is convenient to split 
Hn as follows, 



H 



N+3 



Hn+1 + 



N+l,N+3 



(37) 



for different values of / = 0, 1 and 2 depending on the particular matrix elements one needs 
to evaluate. Hn+i,n+3 contains all the horizontal couplings between the sites + I and 
+ 3 and all the vertical couplings from the rungs + / + 1 until + 3. 
The analogue of eqs.(]3^) are now, 



\H 



f+3 — ^kk, {k — 1, 2, 3) 



Af+3,...,Af+4-fc\<A'fc|-nAr+3-fc,Ar+3|(/^fc/Ar+4-fc,...,Ar- 

N+3{(pl\HN+2,N+3\(p2) N+2,N+3 = ^12\(t>l) N+2 

N+3{4>i\Hn+2,N+3\4>3) N+l,N+2,N+3 = ^13\(t>2) N+l,N+2 + e'13 |0l) Af+1 |0l) Ar+2 

N+3,N+2{(f>2\HN+l,N+3\4>3) N+l,N+2,N+3 = ^23\<Pl) N 



(38) 
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together with their hermitian conjugated. The entries eki can be collected into the symmetric 
3x3 matrix, 



J' 



J + 



J' 



2/2 



V V5^3 



J+J') 



2J 



V 5 



2/2 
5^3 



J+J' 



V 5 



(39) 



while 



-13 



(40) 



The term proportional to e'^^ in (^) is what forces us to introduce the overlap W^q 
The RR's for the energy overlaps are given by, 

EN+i = + eil^Af+2 + '2u{VLi2Em+2,1 + ei2^Af+2,l) 

+2f (^^13-£^Ar+2,2 + ei3-^Ar+2,2 + e'i3W^Ar+2) + ^^^(-EAf+l + ^22Zn+i) 

+v'^{En + eagZAr) + 2M^;(^]23^JV+l,l + e23^iv+i,i) 

-£'Af+3,l = En+2 + eil^7V+2 + u{Vti2EN+2,l + ei2^Ar+2,l) 

+'y(^13-£'iV+2,2 + ei3^Af+2,2 + ^'i?Wn+2 
En+3,2 = ^12-SiV+2,l + ei2^Af+2,l + u{Ej^^i + e222'7v+i) + f (n23-E'Ar+i,i + e23-^JV+l,l) 

The initial data for the energy overlaps are 



(41) 



E, 



0,2 







Eo = Eo,i = 

El = eo J' ^1,1 = eo^' ^^1,2 = (42) 

E2 = 2eo[J' + u{J + J') + u^J] E2,i = eo[2J' + u{J' + J)] ^2,2 = eo( J' + J + 2mJ) 



IV. GENERATING FUNCTION METHODS FOR SOLVING THE RR'S 

The simplest way to find the general term of a series defined by a RR is by introducing 
generating functions. Fan and Ma in have also used generating functions in order to 
find the ground state energy of the dimer state. However in their approach they do not start 
from RR's that generate ground states, and so the appearance of generating functions seems 
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rather obscure. Our method gives a simple and straighforward derivation of their generating 
function methods. Let us illustrate the technique with the derivation of the Binet formula 



A. Number of states 



For this purpose let us define the following generating function, 

F{x) = ^ F^x^ (43) 



N>0 



A simple computation using (p!2|) and (|T^) yields 



Fix) = (44) 

1 — X — 



To recover Fn we can use countour integrals: 

where Co is a countour that encloses the origin counterclockwise. 

For reasons which will become clear later on, it is convenient to perform the change of 
variables x = 1/y, in which case (P5| ) becomes: 

Fn = I ^y^'-'ny) = E Resiy'^-'Fiy)) (46) 

where 

F{y) ^ F(i) (47) 

In (^) Coo is a contour around infinity which picks up all the poles of the function F{y). 
In the example (|4^) we get, 

F{y) = , . (48) 

Noticing that ai,2 = ^(1 i V^) are the two roots of the equation — ?/ — 1 = 0, we get the 
Binet formula ([TT|) 
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We can similarly derive eg . (|22D for the number F^^ of states generated by the RR of 
order u. The generating functions F^'^\x) and F'^iy) are given by: 



2 



X 



F^^Xx) = 1/(1- x-x 

(49) 

F''{y)= y^iy'-y'-'-.-.-y-l) 



B. Norm and Energy of the u = 2 states 

Let us now solve the RR's for Z^v and Ef^f in the u = 2 case. To this end we introduce 
the following generating functions 

(50) 

E{x) = Y.N>0 Enx^, D{x) = Y.N>0 Dnx^ 
From equations (0), (0), (0) and (|0|) we find 

{1 — X — x'^u^)Z{x) = 1 + xuY{x) 

(i-f)r(^)= xz{x) 

(51) 

(1 - X - x^u'^)E{x) = eo{xJ' + 2x'^u'^J)Z{x) + eoXu{2J + J')Y{x) + xuD{x) 
(1 - f )D(x) = x\E{x) + toJ'Z{x) + f (2J + J')tQY{x)] 

Then eliminating Y{x) and D(x) in terms of E{x) and Z{x) we find, 

where A{x), B{x) and C(a;) are polinomials of degrees a = l,b = 3,c = 4 respectively, given 
by 

A{x) = 1 — ^xu 

B{x) = 1 - (1 + - (m^ + lu)x'^ + ^x^u^ (53) 
C{x) = eox [f + xu{2J + J' + 2uJ) - x^u^{\J' + 2uJ) + |xV J] 

We shall show below that eqs.(0) also hold for z/ = 3, where A, B and C are polinomials 
of degrees 3, 6 and 9 respectively. For both, v = 2 and 3, we have the relation. 



a + b = c (54) 
We shall again make the change of variables x = 1/y, defining the new polinomials, 

P{y) = Q{y) = y'B{l/y), R{y) = y''C{l/y) (55) 

In the large limit, and are given by, 

^ Q'ia)' ^ ag'2(c,) 

where Q'{y) = dQ{y)/dy and a denotes the biggest root of the polynomial Q{y). 
The density energy per site of the 2-legged ladder is finally given by, 

n^oo2NZn 2aQ'{a)P{a) ^ ' 



where we have employed eqs.(54j56|). 



C. Norm and Energy of the = 3 states 

A feature of the v = 2 and 3 RR's for the generating functions Zi{x) and Ei{x) , which 
we believe is valid for higher order RR's, is that they can be written in the following compact 
form, 

YX^1qKi^{.x)Z^{x) = e/, 

(/ = 0,...,J/-1) (58) 

where (cq, ei, . . . , Cy^i) = (1, 0, . . . , 0). If Kim{x) is an invertible matrix then the solution of 
these eqs. is given by, 

Zi{x) = Kro\ Ei = Kr^LmnK-^ (59) 

Calling Min the (/,n) minor of the matrix Kin-, then eqs(^9]) have the generic form 
postulated in ( p2D with the identifications. 
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A{x) = Moo, B{x) = detK{x), C{x) = Mok Lm Mi 



10 



(60) 



The matrices K{x) and L{x) for z/ = 3 can be derived from eqs. (|35|) and (^) and they 
read, 



K{x) = 

1 — X — x"^ — v"^ x^ — 2mx (fii2 + 2f Xfi23) 
—X 1 — U xQi2 



—ux 



-X (fii2 + ^23 VX) 



—2vxQis 

V X 

1 



(61) 



iL(x) = 

^ en + ^22xu'^ + x^v{ve^-i + +2uVti2e\.^) 2[uei2 + uve23X + ve[^{x + Qu^isvx'^)] 2vei3 ^ 



(62) 



en + uvx'^Qi2^[^ 
uxe22 



we 12 + 'ye'i3(x + VLi2^i2,vx'^ 
ei2 + vxe23 







V. MINIMIZATION OF THE GROUND STATE ENERGY: RESULTS 

In the previous section we have derived a formula for the ground state energy density, 
eg. (|57[) , in terms of three pohnomials evaluated at the biggest root a of Q{y). The minimiza- 
tion procedure consist now in looking for the minimum of Cqo by varying the parameter u 
for I' = 2, and the parameters u and v for u = 3. In tables 1 and 2 we show the ground state 
energies per site for different values of the coupling constant ratio J/ J', varying through 
strong, intermediate and weak coupling regimes. The values -e^^ / J' and -e'^^ / J' are those 
of the RVA states for v = 2 and 3. The values —e^jJ' are Mean Field values taken from 



, while — e^"/ J' are Lanzcos values taken from 



20 



J/J' 


u 






„Lan / 11 
^oo / "J 








0.375 


0.375000 




0.2 


0.128521 


0.383114 


0.382548 




0.4 


0.323211 


0.40835 


0.405430 




0.6 


0.578928 


0.44853 


0.442424 




0.8 


0.87441 


0.499295 


0.489552 




1 


1.18798 


0.556958 


0.542848 


0.578 


1.25 


1.58519 


0.63518 


0.614473 


0.6687 


1.66 


2.21853 


0.772172 


0.738360 


0.8333 


2.5 


3.39153 


1.06915 


1.002856 


1.18 


5 


5.9777 


1.99285 




2.265 



Tabic 1 



J/J' 


u 


V 






„Lan / 1/ 
^oo 1 











0.375 


0.375000 




0.2 


0.105059 


0.0200031 


0.383195 


0.382548 




0.4 


0.230463 


0.0911438 


0.409442 


0.405430 




0.6 


0.399714 


0.222909 


0.45252 


0.442424 




0.8 


0.626206 


0.42656 


0.507909 


0.489552 




1 


0.91059 


0.716812 


0.571314 


0.542848 


0.578 


1.25 


1.34146 


1.22022 


0.657551 


0.614473 


0.6687 


1.66 


2.19835 


2.41718 


0.808438 


0.738360 


0.8333 


2.5 


4.34045 


6.26652 


1.13384 


1.002856 


1.18 


5 


11.2015 


24.0891 


2.13608 




2.265 



Table 2 



It is also possible to minimize the energy of ladders with finite length, 
show the results obtained with the u — 2> ansatz at the isotropic value J/ J' — 
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In table 3 we 
1. denotes 



the number of vertical rungs. 



N 


u 


V 




3 


0.74225 


0.624464 


0.521564 


4 


0.998703 


0.612503 


0.533225 


5 


0.90358 


0.696107 


0.511372 


6 


0.921508 


0.64184 


0.545944 


7 


0.930216 


0.667329 


0.549671 


8 


0.917952 


0.67691 


0.552377 


9 


0.920783 


0.675904 


0.554453 



Table 3 



Prom these results we can make the following comments 

• For any value of J/ J' there is always a solution which minimizes the energy. In other 
words, the dimer and non dimer RVA states are physical acceptable in the whole range 
of couplings. 

• In the strong coupling regime J/J'< 1 the RVA states give a slightly better ground 
state energy than the mean field result. This later state produces rather unphysical 
results for J/ J' > 1, which does not occur in our case. 

• The nondimer state improves considerably the g.s. energy as compared with the 
dimer state, especially for J/ J' > 1. As J/ J' increase the ratio v/u also increases and 
it is greater than one for J/ J' ~ 1.5. Later on we shall discuss in more detail the 
interpretation of the values of u and v in the isotropic case. 

• In the weak couphng hmit J/ J' » 1 u approaches asymptotically the value 15.94, 
while the energy density approaches - 0.379 J. The later value is much greater than 
the exact result -0.4431 given by the Bethe ansatz solution for the decoupled spin 1/2 
chains. It is however curious that in the limit J/ J' ^ oo the variational parameter u 
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does not go to infinity, yielding a dimer state whose energy is - 0.375 J. This shows 
once more that the resonance mechanism always lowers the energy. 



VI. SPIN CORRELATION LENGTH 

The technique we have used in the previous sections to compute the g.s. energy of the 
RVA states can be easily extended to the evaluation of spin-spin correlators. 



A. Dimer state 

We want to compute the correlator (Si ■ Sj?) between the spin operators at the positions 
1 and R on the ladder ( we assume for simplicity that both sites are on the sublattice A). 
As usual we shall define two auxiliary quantities 



(63) 



gN{R)= (iV|Si-Si^|iV) 
fN{R)= (0i|®(iV-l|Si-Si^|iV) 

From the RR one finds, 

9n+2{R) = gN+i{R) + u^9n{R) + ufN+l{R) 

fN+2iR) = 9N+liR) + ^fN+liR) 

The initial conditions are given by, 

3/m\^ ^ N = R 

g^iR) = 2MR) = <( ^ (g5) 

N <R 



(64) 



These eqs can be derived using the Sutherland's rules |]28| which imply that there is 
only one loop covering giving a contribution to the correlation when the spin operators are 
located at the two boundaries of the chain. Defining the generating functions g{x),f{x) as 
usual we convert the RR's flBSD into, 



[1 — X — x'^v?)g{x) = gR{R)x^ + xuf{x) 

(66) 

(1 - lxu)f{x) = fR{R)x^ + xg{x) 
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Eliminating f{x) we get, 



g{x) = gR{R)x''/B{x) =^ gN{R) ~ a^~''+^gRiR)/Q'{a) (67) 
where a is the biggest root of the polynomial Q{y)- The correlator is finally given for i? >> 1 

by, 

(S.-Snj-^-^.-^m""' (68) 



Zn 2(a-M/2) \2a, 

which gives the expected exponential decaying behaviour with the distance. The correlation 
length is then given by the expression, 

ei/« = — (69) 
u 

From these result we observe that C = exp(l/^) satisfies the cubic equation, 

u^C^ - (2 + u)u^C'^ - (2 + 4:u)u^C + Au^ = (70) 



Setting M = 1 in (|70| ) we obtain the solution of £ = 4.201472 which yields a correlation 
length ^ = 0.696652. This result has been obtained before by White et al. in reference M, 




The correlation length as a function of the ratio J / J' is found by replacing in (|69|) the 
value of u that minimizes the g.s. energy. We collect our results in Table 4. For the isotropic 
case the g.s. corresponds to u = 1.18798, which gives a the correlation length ^ = 0.737253, 
which is slightly bigger than the one given above for u = 1. 



B. i/ = 3 

We shall outline the main steps of the derivation of S,-, since they follow those of the 
dimer case. In particular one gets the same asymptotic behaviour Qn^R) ~ Oi'^~^9R{.R)- 
The main difficulty of the nondimer case is the computation of g^lR). Fortunately we can 
apply some kind of "nested" RR method to compute this quantity. In order to do so we 
need the following definitions, 

gR= {R\Si-Sr\R) 
gR,i= (i?-l|Si-SH|i?) 

gB.,2= (02|® (i?-2|Si-s^|i?) 
^i?„2= (i?-l|Si-S^_i|i?) 

The RR's satisfied by these quantities are given by. 



(71) 



5fij+3 = 2mSi2 fi'R+2,1 + 2t;Si3 {gR+2,2 + gR+2,2) + 2ui;S23 gn+i. 

gR+3,1 = Ul]i2 gR+2,1 + vHvi {gR+2,2 + ^R+2,2) 
gR+3,2 = S12 gR+2,1 + VT,23 gR+1,1 
gR+3,2 = ^12 gR+2,1 + "^^23 gR+1,1 

where we have use the matching eqs. 



(72) 



2(0l|S2|02)l,2 = S12 Si|0i) 



3 



l|S3|03)l,2,3 = Si3 (81 + 82)102)1,2 
3,2(02|83|03)l,2,3 = ^23 8i|0i)i 
3,2(02|82|03)l,2,3 = ^23 Si|0i)i 



(73) 



with 
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1 2 

Si2 = 1/2, Si3 = S23 = -]l - 

The initial data to iterate (|72|) are, 

fl'i = 3/4, 5(11 = 3/4, 5-1,2 = 0, ^1,2 = 

52 = 3m/4, 92,1 = 3n/8, 52,2 = 3/8, 52,2 = 1(1 + ^) 



(74) 



(75) 



Using generating functions we get that the asymptotic behaviour of g^i is given by 
where f3 is the highest root of the following cubic polynomial. 



Q{y) = y'^ - Si2M2/^ - 2f S12 Si3 y - 2t;^Si3 S23 

Hence the correlation length ^ is given by the formula. 



(76) 



(77) 



which contains eq.(^) as a particular case. Indeed setting t> = in ([76|) we get P = u/2. 



In table 4 we show the values of ^ computed from eg. (177|) , for those values of u and v 
that minimizes the g.s. energy. 



J /J' 




^(3) 











0.2 


0.348297 


0.437166 


0.4 


0.471180 


0.608323 


0.6 


0.577543 


0.751286 


0.8 


0.665915 


0.866958 


1 


0.737253 


0.959249 


1.25 


0.807398 


1.04877 


1.66 


0.890758 


1.15205 


2.5 


0.994468 


1.26951 


5 


1.121437 


1.38532 



Table 4 
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Some comments are in order 



There is an improvement of the correlation length specially when J/ J' increases. How- 
ever the value of ^ at the isotropic point is still far from the numerical result which 
is 3.2 1^, [^]. This means that one should consider a RVA ansatz with longer 



bonds. 

If for the dimer state we let u to go to infinity, then a ~ m and we get an upper limit 
for i 

^(2) < i/ln2 = 1.442695 (78) 

Setting M = and taking v large one gets that ~ y2_R/3_ Tj^jg behaviour can be 
understood diagrammatically by constructing the bond configurations that contribute 
to the correlation (Si ■ Sr), which are given by a succesion of states 0i03, and such 
that the loop covering generated by their overlap looks like a braid connecting the two 
extremes of the ladder. The upper bond of ^ is given in this case by, 

^(3) < 3/ln5 = 1.864008, {u = 0) (79) 



VII. DISCUSSION AND PROSPECTS 

The results of these paper give a confirmation of the short range RVB picture of the 
2- legged ladder proposed in 0, especially in the strong and intermediate coupling regimes. 
Indeed, the state 03 can be considered as a pair of two separated topological deffects con- 
nected by a long bond of type /i2 and two dimer bonds in a staggered "high energy" configu- 
ration. The nondimer RVA state generated by this type of local configurations, improves the 
g.s. energy and correlation length of the dimer state, but for J ~ J' one needs to consider 
longer bonds in order to approach the correlation length computed using QMC or DMRG 
methods. 
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It is interesting to evaluate the RVB parameters /ii,/i2 which correspond to the RVA 
parameters u and v. According to eqs.(p!0|, |I8|) they are given, for the case J/ J' = 1, by 

h = 0.954248, h2 = 0.497865 (80) 




FIG. 10. Mixed resonating mechanism involving the building block states of the v = 3 Recur- 
rent Relation. 



The value of /ii ~ 1 is due to the resonance mechanism depicted in figure 3. On the other 
hand the remarkable proximity of /i2 to 1/2 can be interpreted in terms of the resonance 
between horizontal and vertical bonds having NN sites ( see figure 10). This interpretation 
leads us to conjecture that for longer bonds the RVB amplitudes should behave approxi- 
mately for the isotropic ladder as. 

For J = J', and n > 1, /i„ ^ 1/2""^ (81) 

This guess of course agrees with the expected exponential decaying behaviour of bond 
amplitudes of short range states [|^. It would be interesting to confirm (|81|) either by higher 



order RVA ansatzs or by Monte Carlo methods as those of ref p4 . 

The RR method can also be used to compute the value of the spin gap and the string 
order parameter which characterizes the hidden topological LRO of the ladder. These results 
will be presented elsewhere. We have focused in this paper on the two-legged spin ladder 
but the techniques developped so far can in principle be apphed to 4, 6, ... legged ladders, 
and more generally to systems with short range correlations as the spin 1 Heisenberg chains, 
etc. The idea behind the RR method, has some similarities with the RG method of Wilson 



or the DMRG of White p5[ in the sense of constructing the ground state of a system in 
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succesive steps by adding new sites (see also the variational formulation of the DMRG ^^^)- 
The role of u is similar to the role of the number of states m kept in the DMRG method. 
For moderate values of the order u of the RR we can perform analytic computations of the 
g.s. energy and correlations. To improve the accuracy of the RR method would require 
to implement it numerically. In summary the RR method gives us a way to study low 
dimensional quantum systems which is worth to pursue. 
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